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Objectives 

The  objective  of  this  research  continues  to  be  the  study  of  shot-noise  models  and  their  application 
to  the  development  of  computationally  feasible  procedures  for  image- detection  problems.  We  are 
also  interested  in  the  performance  evaluation  of  these  procedures.  These  efforts  are  motivated  by 
applications  to  low-light-level  imaging  as  would  occur  in  low- dose  x-ray  exposures  or  in  night  vision 
systems.  Our  results  are  also  be  applicable  to  photon-limited  optical  communication  systems  and 
to  optical  neural  networks. 

Status  of  Effort 

Considerable  progress  has  been  made  in  understanding  shot  noise.  The  three  major  results  of  our 
research  are  (i)  methods  for  computing  shot-noise  distributions  when  both  the  intensity  function 
and  the  system  impulse  response  are  known;  (ii)  methods  for  estimating  an  unknown  intensity 
when  only  the  impulse  response  is  known;  and  (iii)  methods  for  jointly  estimating  the  intensity  and 
the  impulse  response  when  both  are  unknown  (a  type  of  blind  deconvolution). 

These  results  are  discussed  in  more  detail  below  along  with  other  work  supported  by  this  grant. 

Accomplishments/New  Findings 

•  Shot  noise  that  arises  in  spectroscopy  has  an  underlying  intensity  which  is  a  positive  linear 
combination  of  Voigt  functions.  A  Voigt  function  is  the  convolution  of  a  Gaussian  probability 
density  and  a  Cauchy,  or  Lorentzian,  probability  density.  Because  Voigt  functions  arise 
in  many  different  contexts,  such  as  nuclear  reactor  theory,  atmospheric  transmittance,  and 
spectroscopy,  there  has  been  much  interest  in  computing  them.  In  [1]  (see  the  Publications 
section  below)  a  new  series  for  approximating  Voigt  functions  is  presented,  including  error 
bounds.  An  asymptotic  expansion  for  large  argument  values  is  also  given,  along  with  error 
bounds. 

•  In  many  applications,  one  does  not  know  the  underlying  intensity  of  the  shot  noise,  or  even 
the  system  impulse  response.  For  the  case  when  the  system  impulse  response  is  known,  the 
paper  [5]  uses  a  regularization  technique  involving  a  pair  of  regularization  parameters.  A  new 
method  for  selecting  the  pair  of  parameters  is  presented  and  compared  with  the  minimum 
unbiased  risk  method.  Although  both  methods  recover  the  detail  in  the  intensity  function, 
the  new  method  does  a  better  job  at  suppressing  spurious  oscillations. 

For  the  case  when  both  the  intensity  function  and  the  system  impulse  response  are  unknown 
(blind  estimation),  the  paper  [9]  (which  was  revised  during  95-96)  first  estimates  a  histogram 
of  the  underlying  point  process  and  then  constructs  a  kernel  estimate  of  the  intensity.  The 
estimate  of  the  unknown  impulse  response  is  constructed  via  a  regularized  back-substitution 
of  a  discrete-time  convolution  with  the  estimated  histogram. 

•  Even  if  the  shot-noise  intensity  and  system  impulse  response  are  known,  it  is  quite  difficult  to 
compute  the  cumulative  distribution  and  density  for  shot  noise.  Being  able  to  compute  these 
functions  will  enable  the  construction  of  detectors  and  the  evaluation  of  their  performance. 
In  [6],  the  computation  of  the  cumulative  distribution  (cdf),  the  complementary  cdf  (ccdf), 
and  the  density  of  certain  shot-noise  random  variables  is  discussed.  After  subtracting  off 
a  few  terms  that  can  be  computed  in  closed  form,  what  remains  can  be  approximated  by 
a  general  method  for  approximating  samples  of  a  cdf  or  ccdf  by  summing  a  Fourier  series 
whose  coefficients  are  modulated  samples  of  their  characteristic  function.  To  approximate 
the  density,  a  spline  is  fit  to  the  cdf  samples  and  then  differentiated.  When  the  density  has 
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corners,  it  is  important  that  the  spline  have  coincident  knots  at  these  locations.  For  shot-noise 
densities,  these  locations  are  easily  identified. 

In  95-96,  we  have  studied  Monte-Carlo  estimation  of  shot-noise  cdfs  for  comparison  with  the 
above  method.  Although  Monte-Carlo  methods  require  more  computation,  we  found  that 
they  could  be  speeded  up  using  a  change  of  measure  in  which  the  underlying  Poisson  process 
is  simulated  with  an  alternative  intensity. 

The  methods  of  paper  [6]  suggested  a  similar  approach  could  be  used  to  recover  the  cdf  and 
ccdf  of  a  counting  distribution  from  its  characteristic  function.  This  result  is  reported  in  [8] 
and  was  used  to  evaluate  the  shot-noise  limited  performance  of  an  optical  neural  network  [7]. 

•  Shot-noise  with  random  impulse  response  arises  in  the  study  of  avalanche  photodiodes.  The 
effects  of  dead  space  on  the  bit-error  rate  for  an  optical  receiver  have  been  studied  in  [2].  In 
particular,  it  is  found  that  dead  space  enhances  the  performance  at  relatively  low  data  rates, 
and  degrades  performance  at  high  data  rates.  Thus,  in  order  to  correctly  predict  performance 
of  high-data-rate  systems,  it  is  critical  to  account  for  the  dead-space  effect. 

•  Other  research  supported  by  this  grant  included  an  investigation  into  wavelets  for  discrete¬ 
time  periodic  signals  [4],  and  the  capacity  region  of  the  multiple- access  arbitrarily  varying 
channel  [3]. 
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